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ABSTRACT 

Galaxy pairs provide a potentially powerful means of studying triggered star formation from galaxy interac- 
tions. We use a large cosmological N-body simulation coupled with a well-tested semi-analytic substructure 
model to demonstrate that the majority of galaxies in close pairs reside within cluster or group-size halos and 
therefore represent a biased population, poorly suited for direct comparison to "field" galaxies. Thus, the fre- 
quent observation that some types of galaxies in pairs have redder colors than "field" galaxies is primarily a 
selection effect. We use our simulations to devise a means to select galaxy pairs that are isolated in their dark 
matter halos with respect to other massive subhalos (N= 2 halos) and to select a control sample of isolated 
galaxies (N= 1 halos) for comparison. We then apply these selection criteria to a volume-limited subset of 
the 2dF Galaxy Redshift Survey with Mb.j < -19 and obtain the first clean measure of the typical fraction of 
galaxies affected by triggered star formation and the average elevation in the star formation rate. We find that 
24% (30.5 %) of these L* and sub-L* galaxies in isolated 50 (30) h"' kpc pairs exhibit star formation that is 
boosted by a factor of > 5 above their average past value, while only 10% of isolated galaxies in the control 
sample show this level of enhancement. Thus, 14% (20 %) of the galaxies in these close pairs show clear trig- 
gered star formation. Our orbit models suggest that 12% (16%) of 50 (30) h"' kpc close pairs that are isolated 
according to our definition have had a close (< 30 h"' kpc) pass within the last Gyr Thus, the data are broadly 
consistent with a scenario in which most or all close passes of isolated pairs result in triggered star formation. 
The isolation criteria we develop provide a means to constrain star formation and feedback prescriptions in hy- 
drodynamic simulations and a very general method of understanding the importance of triggered star formation 
in a cosmological context. 

Subject headings: cosmology: theory, large-scale structure of universe — galaxies: formation, evolution, high- 
redshift, interactions, statistics 



1. INTRODUCTION 

Galaxy interactions and mergers help drive galaxy evo- 
lution. In the concordance ACDM cosmological model, 
nearly every galaxy has had a major merger at least once 
over cosmic time (e.g., iMaller et alj l2006t IStewart et alj 
l2007h . Major mergers and interactions consume available 
gas, producing ste llar populations wh i ch sub s equently redden 
with age (e.g., iLarson & TinslevI Il978t iKennicutt et alj 
[T987h . Major mergers destroy disks, tur ning galaxies into 
spheroids (e.g., iToomre & Toomrj Il972h . but also poten- 
tially cre ating large disk galaxies in gas-rich scenarios 
(e.g., Rob ertson et al.l 1200 6). Close galaxy passes and 
minor mergers are even more frequent. These perturba- 
tions send gas into the centers of galaxies, where they 
may contribute to bulge comp o nents a nd even feed black 
holes (e.g., Kennicutt & Keel 1984; Barnes & Hernquist 
1992; Mihos & Hernquist 1996; Barton et al. 2000; Combes, 
2001; Barton etal. 2001; Barton Gillespie et al. 2003t 
Kannappan, Jansen, & Barton 2004; Kormendv & Kennicutfl 
^04; Freedman Woods, Geller, & Barton 2006; Lin et al. 
2007t lAlonso et al.ll2007t iFreedman Woods & Gelleil 120071) . 
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In the hierarchical picture of galaxy formation, mergers 
play an even more important role during earlier epochs 
than they do t oday, possibly tr i ggerin g the high star for- 
mation rates flLowenthal et alj [l997t iKolatt et al] 119991; 
ISomerville eTal] 120011; IWechsler et al] 1200 Ih observed in 
some "Lyman br eak" galaxies at z > 3 dSteidel & Hamiltonl 
ll992l:[Steidel et a l. 1996) or luminous submillimeter galaxies 
(e.g..' Chapmanet al.ll2005l) . 

Existing detailed studies of galaxy interactions show that 
triggered star formation and morphological evolution can be 
very rapid and intense. However, we can only observe a 
"snapshot" of the evolution of galaxy pairs. With incom- 
plete phase-space information, we do not know their true fre- 
quency, evolutionary timescales, orbits, or fates. Thus, we 
have neither a complete understanding of the frequency of 
interactions and mergers nor knowledge of their impact on 
galaxy evolution. Many processes such as ram pressure strip- 
ping, interaction with a cluster or group potential, or high- 
speed "harassment," may also be very important in estab- 
lishing the morphologies and stellar populations of galaxies, 
and in driving the relationships between morphology, color, 
and local environmental density (e.g., Gunn & Gott 197 J; 
Dressier '1980; Post man & Gelleii 119841: iMoore et all 119961; 
Blanton et al. 2005a). 

Hydrodynamic simulations of galaxy interactions can make 
predictions for the expected star formation rate from any 
type of galactic colhsion (e.g., [B arnes & Hern quist! 119921; 
Mihos & Hernquist 1996"; Springel, Di Matt eo. & Hernquisj 
2005; Cox et al. 2006; Perez et al. 2006a b). However, these 
predictions are sensitive to prescriptions for baryonic physics 
that are highly uncertain. In particular, hydrodynamic simu- 
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lations must rely on effective models for gas dynamics, cool- 
ing, and star formation that these simulations cannot hope to 
resolve. Uncertainties in these effective theories are an enor- 
mous problem for galaxy formation theory and detailed tests 
of these prescriptions are required. Unfortunately such tests 
are not simple; the full phase-space information is not known 
for many real systems. However, empirical, statistical studies 
of galaxy interactions and mergers — when combined with 
predictions of galaxy orbits in a now well-established cosmo- 
logical model — hold promise for uncovering the role of in- 
teractions and mergers in galaxy evolution. 

The first step in studying the effects of interactions and 
mergers is to identify systems that are undergoing these pro- 
cesses. Using morphological distortion as an indicator of 
galaxy interactions does yield a subset of the pairs that have 
definitely had a close pass. However, relying on tidal distor- 
tion leads to missed interacting pairs because (1) sensitivity 
to low-surface-brightness features is a strong function of see- 
ing, depth, and redshift, (2) the morphological features are 
short-lived (^100 Myr), and (3) many of the features, like 
tidal tails, are resonanc e effects th at only appear in prograde 
disk galaxy encounters jToomre & T oomre 1972). For a more 
complete census of the effects of interactions, it is important 
to select pairs or systems of galaxies based only on proximity 
in redshift and projected separation, then to use mock catalogs 
constructed from cosmological simulations to understand this 
selection. 

Historically, studies of distorted galaxies and galaxies in 
pairs have been very revealing. The early stages of galaxy 
interactions can drive gas into the center of a progenitor 
galaxy and trigger an early episode of central star formation 
long before the final merger (e.g ., Larson & Tinslev 1978; 
Joserih et all 1984 iKeiinicutt et al.lll987: Mihos & Hernq uist. 
1994 Il996l) . The strength of the optical emission line as- 
sociated with this star formation correlates with the sepa- 
ration jaf the pair on the sky (AD) and in redshift (AV; 
fearton et ai.ll2000l) . The optically-detected star formation is 
strongest in the central few kpc, often dominating the opti- 
cal light; in the case of a late-type spiral, it c an also con- 
tribute substantially to the formation of a bulge ( Tissera et al.l 
12002; Barton Gillespie et al. 2003). This process may be 
the primary mechanism for the forma ti on of late-type bulges 
(iKannappan. Jansen. & Bartonll2004h . iKewley et al] (l2006h 
show that the metalUcities of these galaxies provide a "smok- 
ing gun" for gas infall from the outskirts of the disks. Galax- 
ies with (optically) strong central starbursts have metallicities 
that are lower than average for their luminosities, consistent 
with a starburst occurring in gas that was driven into the nu- 
cleus from the metal poor outskirts of the progenitor's disk. 

New, large red shift surveys such as the 2dFGRS 
(IColless et alJ2001h and Sloan Digital Sky Survey (lYork et all 
12000') provide a means of exploring ever-larger samples of 
galaxies in pairs, although these samples must be approached 
with caution because mechanical spectrograph constraints 
make them deficient in close pairs. The correlations between 
orbital parameters and star-forming properties of galaxies in 
pairs have been verified in both the 2dFGRS (Lambas et ajj 
l2003h . and the SDSS (N ikolic et al. 2004; Luoetal. 200d). 
However, the studies also appear to reveal unexpected phe- 
nomena including apparently suppressed star formation in 
widely-separated pai rs and red tails in the color distributions 
of paired galaxies ( Lambas et al.l l2003t lAlonso et al.l l2004l 
12006; Luo et al. 2007^^ 

Using these empirical results to arrive at a true measure of 



the amount and timescales of triggered star formation is dif- 
ficult. Pair samples are complex: they include some inter- 
loper s and are ri ch in galaxies i n a va riety of environments 
(e.g.. lAlonso et al. 2006; SoaresI 120071) . The progenitors of 
galaxy inter actions come from a mixture of galaxy types (e.g., 
iFocardi etal. 2006). In dense environments, many will have 
already experienced multiple close passes and mergers. They 
may have consumed or lost much of their gas, leaving little to 
form new stars. However, in sparser environments, the pro- 
genitors of the interaction may have had very few previous 
interactions. They may have large remaining gas reservoirs. 
Thus, the efficacy of an interaction in triggering star formation 
should depend in detail on the environment of the interacting 
galaxies. 

Dense clusters of galaxies are straightforward to identify 
with confidence. However, the low-speed encounters that trig- 
ger star formation are extremely unlikely in these clusters. An 
understanding of tidally-triggered star formation requires ac- 
curate probes of isolated galaxies and sparse loose groups. 
Unfortunately, accurate environmental statistics are difficult 
to interpret precisely in this regime. Group-finding algorithms 
are ef fective at investigating environments ( Huc hra & Gelled 
119821) . However, even highly tuned group-finding algo- 
rithms construct f alse groups and rn i ss group members (e.g. , 
'Yang et al."2005'; Gerke et al. 2005'; Weinm ann et"ani2006t 
£erlind et al. 2006; Koester et al. 2007). Nearest-neighbor 
and counts-in-cylin der statistics also have a large scatter (e.g., 
iBerrier et al.1l2007h . 

The use of mock galaxy catalogs based on numerical sim- 
ulations of structure formation in the standard cosmology can 
help sort these issues out. Direct numerical simulation is dif- 
ficult because obtaining the necessary resolution to be com- 
plete in close pairs while simultaneously modeling a large 
enough volume to reduce sample variance is computationally 
expensive. An alternative and proven method is to use an ana- 
lytic model for dark matter halo substructure (so-called "sub- 
halos". e.g.. IZentner & Bullock! 120031: iTaylor & BabulH2004t 
IZentner et al.l l2005b in conjunction with an A^-body simula- 
tion of a cosmologically -relevant volume. The analytic model 
treats halo substructure using a method that is approximate, 
but free of inherent resolution limits, and extends the effective 
numerical resolution in the simulated volume, ensuring com- 
pleteness in the dense environments where many close pairs 
reside. The analytic model also allows for a quantification of 
shot-noise contributions to close-pair samples that may arise, 
for exampl e, from the particula rs of galaxy orbits and may be 
significant (iBerrier et al.l l2006). Such models have been de- 
veloped and validated by comparison to direct A^-body simu- 
lations i n the regimes where the two techniques are commen- 
surable (IZentner et al.l2005 h. Thus, the time is ripe to develop 
new methods aimed at understanding the detailed make-up of 
pairs of galaxies selected from redshift surveys. 

We study mock catalogs constructed using the hybrid ap- 
proach desc ribed briefly i n the p revio us paragraph and ex - 
tensively in IZentner et all (l2005l) and IBerrier et all (l2006l) . 
This approach has already succeeded in explaining the ap- 
parent lack of evolution in the close pa ir fraction ob served 
in inte rmediate-redshift surveys (iLin et al . 2004; Berrier et al.l 
|2006|) . In the present study, we focus on interpreting the star- 
forming properties of galaxies in pairs. With the complexities 
of the large-scale environments of pairs in mind, we explore 
the construction of appropriate samples and control samples 
for galaxies in pairs. In particular, we examine ways to isolate 
the immediate effects of triggered star formation from other 
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environmental processes. Here, we apply the analysis to a 
volume-limited sample of galaxies in the 2dFGRS; however, 
the techniques we discuss are generally applicable to other 
studies of tight sub-groupings of galaxies. 

In §|2]we describe the numerical methods and the observa- 
tional data that we use. §|3]contains a description of the model 
predictions for the environments of galaxies in pairs; in this 
section we examine the problem of trying to construct a con- 
trol sample of objects that are not in pairs. In §|4l we restrict to 
the most isolated pairs and individual galaxies and show that 
isolated pairs are almost purely dark matter halos containing 
N = 2 galaxies inside their virial radii. The isolated galaxies 
we select are an appropriate control sample for the progeni- 
tors of an interaction. We examine galaxies in pairs in t he 2dF 
Galaxy Redshift Survey (2dFGRS Coll ess et al.ll200TI) in § |5] 
and measure the amount of triggered star formation in isolated 
pairs relative to the control. We §[7] contains a brief descrip- 
tion of the situation in more complex environments, and we 
conclude in § [8] In future papers, we plan to explore more 
complex environments in greater detail. The ultimate goal of 
this and similar studies is to isolate the effects of interactions 
and then re-integrate the results into a complete cosmological 
interaction history for galaxies, thus measuring the amount of 
galaxy evolution triggered by interactions and mergers. 

2. THE MODELS AND THE DATA 

The first step in our study is to construct mock catalogs of 
galaxies and study the environments of galaxies in projected, 
close-pair configurations in order to understand the selection 
of such objects. We use the information gleaned from the 
mock samples to aid in the interpretation of close pairs in the 
2dFGRS. In this section, we describe the models used to pro- 
duce mock catalogs and the observational data in turn. 

2.1. Mock Galaxy Catalogs 

Our model begins with a cosmological N-body simulation 
from which we extract positions and masses of host dark mat- 
ter halos within a cosmological volume. We populate the hos t 
halos with subhalos using the method of lZentner et al.l (l2005h . 
Our method for producing mock catalogs has been discussed 
in detail in Berrier et al. (2006 ). We provide a br ief overview 
here, and refer the reader to iBerrier et al I (l2006l) for a more 
detailed discussion. 

Our numerical simulation was run with the Ad aptive Re- 
finement Tree N-body code (" Kravtsov et al.l [19971) in a stan- 
dard ACDM cosmology with Q,„ = 0.3, Qa = 0.7, Hq = 70 km 
s"' Mpc~', and as = 0.9. The simulation followed the evolu- 
tion of 512^' particles in comoving box of 120 h"' Mpc on a 
side, implying a particle mass ofm^ 2± 1.1 x lO^h"'M0. The 
simulation grid was refined down to a minimum cell size of 
/ipeak — 1-8 h~' kpc on a side. The simulation was previously 
discussed i n Tasitsiomi et al. (2004), Zentner et al.l (l2005l) . 
lAIlpood et al.. (.2006,) . and ,Wechsler et al. (2006.) . Host ha- 
los in the simulation were identified usi ng a variant of the 
Bound Den sity Maxima algorithm (BDM.'Klypin et al.l l 19991: 
lKravtsov~et al. 2004). We define a halo virial radius as the ra- 
dius of the sphere, centered on the density peak, within which 
the mean density is A,.(z) times the mean density of the uni- 
verse, pm- The virial overdensity Ay(z), is given by the spher- 
ical top-hat collapse approximation and we com pute it using 
the fitting function of Bryan & Norman] (Il998l) . Host halos 
are identified as those halos whose centers do not lie within 
the virial radius of another halo. 

In order to determine the substructure properties in each 



host dark matter halo, we model their mass accretion his- 
tories and track their substructure content using an analytic 
technique that exploits numerous simplifying approximations 
and several sc aling relations derived from direct simulation 
(IZentner et a l. 2005). For each host halo of mass M at red- 
shift z in our simulation volume, we rand omly generate a mass 
accretion history using the method of ISomerville & Kolata 
(Il999l) . At each merger event, we assign an initial orbital en- 
ergy and angular momentum to the inf ailing object accord- 
ing to the probability distributions for these quantities de- 
rived from cosmological A^-body simulations (IZentner et al.l 
l2005h . Each accreted system becomes a subhalo at the 
time it is accreted. It is assigned a mass and a corre- 
sponding maximum circular velocity at this time, Vinax= Vin 
. We integrate the orbit of the subhalo in the potential of 
the main halo from the time of accretion until z = 0. We 
model both tidal mass loss and internal heating of subha- 
los as well as dynamical friction using a modified form of 
the Chandrasekhar formul a ( Chandrasekhar 1943) suggested 
by [Hashimoto. Funato. & Making (.2003). As subhalos orbit 
within their hosts, they lose mass and their maximum circular 
velocities decrease as the profiles are heated by interactions. 
We remove galaxy subhalos from our catalogs once their max- 
imum circular velocities drop below Vmax= 80 km s~' . This 
rough criterion is used to mimic the dissolution of the observ- 
able galaxy as a result of these interactions. 

This procedure produces a population of subhalos within 
each host dark matter halo in the volume. In fact, the proce- 
dure is statistical and relies on realizations of the small-scale 
density field and the orbital parameters for infalling struc- 
ture. As such, each host may be assigned numerous different 
subhalo populations that differ in their detail due to realizing 
these statistical distributions with finite samples. We include 
such variations by producing four mock catalogs from four 
independently-realized subhalo populations for each host. 

The next step is to map galaxies onto the host halo and sub- 
halo populations in the model. Each subhalo has an associ- 
ated Vin circular velocity that acts as a proxy for its luminos- 
ity and each host halo has a corresponding "central galaxy" 
with a Vin= Vmax- We assign luminosities to halos by match- 
ing volume number densities of galaxies between the simula- 
tion and the data; thus, we ass ume that luminosity is a mono- 
tonic function of In Berrier et al.l (12006.) . we compared 
the close pair fraction predicted by this method to close pair 
counts in the UZC and found good agreement if we associ- 
ated galaxy luminosity in a one-to-one way with Vin. We use 
th is association throughout this paper. In addition to the work 
of lBerrier et al.l (i200 6) which shows the two-point correlation 
function extending down to ^^50 h"' kpc, a similar model has 
also been shown to reproduce the larger scale galaxy two- 
point correla tion function as a function of luminosity, scale, 
and redshift (Conroy et al. 2006). 

As discussed below, we focus on halos with Vin > 160.5 
km s~' in order to define a mock galaxy catalog with the 
same number density as the Mbj < -19 population we con- 
sider from the 2dF survey. The numerical model predicts the 
number of galaxies with Vin > 160.5 km s"' that occupy ev- 
ery host halo in the simulation, N. As discussed below, most 
galaxies in the simulation reside in halos by themselves, with 
N = 1, but galaxies in highly-clustered regions tend to reside 
in massive host halos with multiple galaxies, N > 9. In what 
follows we analyze the full simulation volume; thus, all re- 
sults are appropriate for comparison to a volume-limited sam- 
ple. In addition, we use all four mock catalogs constructed 
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from the four, independent realizations of halo substructure 
for each host unless otherwise stated. 

2.2. Pairs in the 2dF Survey 

We examine star formation in pairs se lected from the public 
2dFGRS database (iCoUess et alj|2001h. In contrast to p revi- 
ous studies of 2dFGRS pairs (e.g.. IXlonso et al.]|2006h . we 
construct a volume-limited sample with galaxies to Mbj = - 
19, assuming 17a = 0.7, n„, = 0.3, and Hq = 70 km s"' Mpc"'. 
The redshift range is from 0.010 to 0.0875. We restrict the 
study to two simple rectangles in the sky, with coordinates 
148.97" < a < 209.42", -4.78" <5< 2.21" and 29.68" < a 
< 50.94", -34.18" <S < -25.32". While these strips are rel- 
atively complete, there are some omitted regions near bright 
stars. We only select pairs from the portion of the sample not 
within 700 h~' kpc of the edges of the rectangles on the sky 
or within 1000 km s"' of our redshift limits in order to probe 
the full environments of our targeted objects. 

The full sample covers a solid angle of 0.317 steradians and 
a volume of 1.79 x lO*" h~^ Mpc^. Integrating the luminos- 
ity function reveals an expected object density of 0.0107 ± 
0.0008 galaxies h^ Mpc'^ to Mbj = -19 (ICroton et al.ll2005h : 
this density corresponds to a cutoff halo circular speed Vin = 
160.5 km s"' in the simulation, which we use consistently 
throughout the paper. The targeted sample excluding the 
edges includes 22,601 galaxies, with 1,344 galaxies in close 
(< 50 h-i kpc) pairs. 

The primary problem faced in our data analysis is survey 
incomplet eness. Overall, t he 2dFGRS is approximately 86% 
complete ( Cross et al]|200 4). In a single observation, the 2dF- 
GRS cannot distinguish objects closer than 25 arcseconds 
on the sky. At the minimum, median, and maximum red- 
shift of galaxies in the volume-limited sample we consider, 
this fiber collision separation corresponds to 4.8, 23, and 28 
h~' kpc, respectively. However, repeated measurements of 
the same field with different fiber confi gurations allow red- 
shift measurement s for some close pairs (iLambas et ani2003l; 
lAlonso et al.ll2004l) . In dense clusters, where galaxy pairs are 
preferentially found, there are often many more objects than 
fibers, leading to a differential incompleteness with environ- 
ment that has been c haracterized by the 2dFGRS survey team 
(IColless etal.ll2001h . We use their tools to probe the effects of 
incompleteness in §|5] 

Although it is not a primary focus of this work, we also 
verify the results of our st udy using the Sloan Digital Sky Sur- 
vey DR4 ( Adelman-McCarthv et al. 2006) NYU Value-Added 
Galaxy Catalog data (Blanton et al. 2005b). We construct 
volume-limited catalogs to Mr = -19 + 51og/i and Mr = -20 + 
Slogh and consider both galaxies with measured and unmea- 
sured redshifts. We use star formation indicator s including 
the star formation rate per unit mass from Brinc hmann et al.l 
(l2004l) . In general, the SDSS results are in qualitative agree- 
ment with the 2dF results we focus on here, but the use of 
r-band selection and of different star formation indicators in- 
troduce quantitative differences. 

3. THE ENVIRONMENTS OF GALAXY PAIRS 

To isolate triggered star formation from star formation that 
is not triggered by interactions, we must compare a well- 
defined pair sample to the appropriate control sample — typ- 
ical examples of the immediate progenitors of the interaction. 
What types of galaxies appear in pairs? To answer this ques- 
tion, we investigate the environments of pairs selected in our 
simulated volume-limited redshift survey. 
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Fig. I. — The fraction of galaxies in the simulation with a total of N 
galaxies in their host halos as a function of Dn, the projected distance to its 
nearest neighbor. For each galaxy in the mock catalog, we measure Dn within 
Ay= 1000 km s"', and the total number of galaxies in its host halo, N. The 
different lines correspond to different bins of N, N = 1 exactly (solid black), 
2 (dotted green), 3 (dashed magenta), 4-8 (long-dashed cyan), and > 9 (dot- 
dashed blue) galaxies in the host halo. We shade the "pair zone," galaxies 
with neighbors within 50 h"' kpc. Close pairs are preferentially found in 
very well-populated halos. 

For each simulated "galaxy" (halo or subhalo) above our 
cutoff circular velocity in the model, we compute Dn, the 
projected distance to the object's nearest neighbor with Ay< 
1000 km s"' . We also measure N, and the total number of sim- 
ulated galaxies that lie within the same host dark matter halo 
as the object. Typically, close pairs are defined based on small 
projected separations, A D < 30- 100 h~' kpc , with 50 h~' kpc 
a commonly used value (iBarton et al.ll200(it iLin et al.ll2004t 
iBerrier et all2006l) . within AV= 1000 km s"' of its neighbor. 
Because we have full three-dimensional information from the 
simulation, we can explore the true nature of the "apparent" 
pairs selected with this technique. 

We begin our exploration of close pair selection in Figure[T] 
To construct this figure, we first compute Dn and N for each 
simulated galaxy, using the periodic boundary conditions, if 
necessary, to fully sample its environment. We then split the 
objects in each bin of Dn into subsets based on the multiplic- 
ities of the halos in which the objects reside, A'^a/(N|DN). The 
ranges of the N bins are shown in the labels of Fig. [T] Fi- 
nally, at each value of projected nearest-neighbor separation, 
we computed the fraction of galaxies within that bin that are 
in the different subsets of N. Fig.[T]shows this fraction among 
five different bins of N as a function of Dn. In each Dn bin, 
the sum of all fractions is 1 . For reference, the shaded region 
indicates the 50 h"' kpc "pair zone." We note that the x-axis, 
Dn, is the directly measurable quantity from redshift survey 
data, while N is not directly measurable. 

Fig. [T] immediately reveals a major issue in selecting ap- 
propriate control samples for galaxies in pairs. Pairs in the 
simulation reside in a highly skewed distribution of environ- 
ments. In particular, pairs are vastly over-weighted toward 
galaxies that reside within well-populated host dark matter 
halos compared to the average. For example, galaxies with 
N > 9 comprise 39% of the pair sample but only 19% of 
the simulation as a whole. In contrast, while isolated (N=l) 
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Fig. 2. — The normalized distribution, dn/dlog(M), of host halo masses 
for 50 h"' kpc pairs (dotted red), closer 30 h"' kpc pairs (dashed blue) and 
the full sample (solid black) in the simulation. Pairs reside in much more 
massive halos than typical "field" galaxies. 

galaxies make up 56% of the simulation as a whole, only 
3% of the apparent close pairs are actually isolated galaxies. 
These isolated galaxies i n apparent pa i rs are interlopers ac- 
cording to the analysis of iBerrier et al.l (l2006t) . The abundant 
pairs in cluster and large-group systems are not necessarily 
imminent mergers, nor ar e they nec essarily even interacting 
directly with each other iBailin et al. (2007) note a similar 
problem in the study of galactic satellite populations. 

From a qualitative perspective, Figure[T]explains many past 
observations regarding the star-forming properties of galax- 
ies in pairs. For example, pair samples previously analyzed 
in the 2dFGRS show an increased star formation rate rela- 
tive to control samples at close separations (AD< 30 h ~' kpc) 
where the interaction has a dominant effect (.Lambas et al.l 
l2003h . However, at larger separations where the interaction 
is less important (100 < AD < 200 h"' kpc), the typical pairs 
exhib it less star formatio n than the field (e.g., Lambas et al. 
120031: lAlonsoet aLll2004 l2006l) . In our model, these more 
widely-separated pairs are also in denser environments than 
the field: 36% are in N > 9-galaxy halos while only 19% are 
isolated. This depressed star formation in widely-separated 
pairs relative to the field results from the fact that galaxies liv- 
ing in more massive and populated systems have suppressed 
star formation on average. 

Because pairs are preferentially found in denser environ- 
ments, the naive comparison of star formation rates between 
pairs and typical field galaxies will artificially underestimate 
any elevation in star formation rate that is directly triggered 
by an interaction. Fig. |2] shows the distribution of host halo 
masses for galaxies in close (AD< 50 and 30 h"' kpc) pairs 
and the sample as a whole. This Figure emphasizes how 
skewed the environments actually are. The average host halo 
mass for the full sample is < Mhost >= 3.2 x lO'"* h"' Mq; 
for the close (50 h"' kpc) and closer (30 h"' kpc) samples 
of pairs it is 6.0 and 5.7x10'^ h"' M0, respectively. More- 
over, the mode of the full sample, ~ lO'^ h"' Mq, is almost 
completely depleted in the pair sa mple. The analysis of the 
SDSS by WeTn mann et al. I (|2006) tracks the dependence of 
galaxy properties on host halo mass. We can use their results 
to understand the implications for galaxies in pairs, which 
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Fig. 3. — The distribution of host halo mass as a function of distance 
to the nearest neighbor, Dn, in the simulation for a single realization of the 
substructure model. (Left) For galaxies with at least one luminous neighbor 
within 700 h"' kpc and AV= 1000 km s"', we plot the log of the host halo 
mass as a function of projected distance to the nearest neighbor in the simula- 
tion. The point colors and types segregate galaxies based on halo occupation. 
We plot single-galaxy halos (solid black triangles), N = 2 (open green cir- 
cles), 3 (filled magenta circles), 4-8 (open cyan triangles), and > 9-galaxy 
halos (blue crosses). The solid red line is the virial radius of halos; the dotted 
red line is the virial radius reduced by a factor of ^2/3 to account for the 
typical difference between projected and actual separation. (Right) We also 
plot the host mass histogram of the host halo of simulated galaxies with no 
neighbors within 700 h"' kpc. 

are preferentially selected from higher-mass halos. Accord- 
ing to their results the "average" color of a galaxy in an "av- 
erage" halo hosting a pair would be ^0.05 magnitudes red- 
der in g-r color, have a specific star formation rate (e.g., 
iBrinchmarm et al. 2004) that is ^20% lower, and have a late- 
type fraction that is ^5% lower than the corresponding av- 
erages for the field. Thus, the suppressed star formation ob- 
served in some widely-separated galaxy pairs is not related to 
a single interaction, but is a selection effect associated with 
cluster and group processes that suppress star formation. This 
selection effect must be accounted for in a fair analysis of 
triggered star formation. 

Fig.[3]illustrates the nature of the differences between halos 
that host close pairs and "typical" dark matter halos. This 
figure plots galaxy host halo mass as a function of projected 
separation to its nearest neighbor up to 700 h"' kpc and within 
1000 km s"' in our mock catalogs. The colors and point types 
segregate host halos based on N, the number of galaxies they 
host. The structure of this plot provides a very general means 
of viewing nearest-neighbor statistics. The isolated halos line 
up along the lower part of the plot; their locus is analogous to 
the "two-halo" term in the correlation function insofar as the 
pairs are projections of two coincidentally aligned host halos. 
The more massive and populated dark matter halos line up 
along the left vertical edge of the plot. Again, pairs are an 
extremely skewed population. In contrast, galaxies that are 
isolated within ^400-500 h"' kpc are almost purely the only 
occupants in their lower-mass halos. 

In summary, a comparison of the properties of typical pairs 
and typical isolated galaxies selected from the universe as 
a whole is not the appropriate comparison for isolating star 
formation triggered by galaxy interactions. The "field" is 
dominated by isolated, low-mass host halos with a mix of 
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a few more luminous systems. Pairs are skewed, preferen- 
tially residing in higher-mass host halos. Thus, comparisons 
of the star-forming properties of pairs and the "field" do not 
isolate the effects of the recent interaction. The properties 
of close pairs reflect all of the other processes that occur in 
dense group and cluster environments that may act to sup- 
press star formation during late-time interactions. As such, 
naive comparisons will dramatically underestimate the rise in 
star formation rate triggered by galaxy interactions. In the 
next section, we describe one means of constructing samples 
that better isolate the effects of galaxy interactions from other 
environment-related processes. 

4. ISOLATING TRIGGERED STAR EORMATION 

Many different measures of large-scale environment have 
been developed to understand the properties of galaxies as a 
function of their surroundings. With the advent of redshift 
surveys, these studies ha ve fo cused on group-finding algo- 
rithms (Huchra & Geller 19821) . «th nearest neighbor statis- 
tics (Dressier 1980), and local galaxy count measures (e.g., 
iBlanton et al. 2005a; Berrier et al. 2007). Here, we adopt 
N700, the number of neighboring galaxies within 700 h~' kpc 
on the sky and AV= 1000 km s~' in redshift of the galaxy 
in question (see lBerrier et al.ir2007h . We note that our results 
are qualitatively insensitive to choices of environment statistic 
scale in the range of 700-1000 h"' kpc, but 700 h"' kpc yields 
a large enough sample of galaxies in pairs in the 2dFGRS 
(§|5]l. The use of other statistics such as group-finding algo- 
rithms would yield somewhat different results, but the basic 
problem of separating dense systems from nearby but unasso- 
ciated isolated dark matter halos remains. 

Choosing galaxies with limited numbers of companions on 
~1 Mpc scale environments is an effective way to preferen- 
tiall y select galaxies in different types of dark matter halos 
(see lFocardi et al.|[200 6'). Fig.|4]shows the distribution of pro- 
jected Dn as a function of host halo mass for the the restricted 
set of environments corresponding to galaxies with N700 = 1 . 
These are galaxies that have exactly one companion within 
700 h"' kpc and 1000 km s"'. This restriction is extremely 
clean and effective. Widely-separated pairs are almost exclu- 
sively isolated halos and the close pairs are almost purely N 
= 2 systems. Fig. |5] shows the fraction of host halos of each 
type as a function of separation in this restricted sample and 
the corresponding host halo masses in the pairs and full dis- 
tributions. The fraction of N = 2 hosts goes from over 90% 
at the closest separations to a 50/50 split of N = 1 and N = 2 
halos at ^175 h~' kpc, to nearly 100% N = 1 halos beyond 
400 h"' kpc. 

The key to the use of ultra-low-density environments is that 
the restriction provides an extremely clean sample of pairs in 
N = 2 halos at small separations and a corresponding con- 
trol sample of isolated galaxies. This allows us to remove 
any effects that may be due to numerous and repeated inter- 
actions in dense cluster environments and study the relative 
effect of individual interactions. The sample average masses 
are < Mhost >=4.2 and 3.9 x lO'- h"' Mq, respectively, for 
the objects that could be identified as isolated 50 h"' kpc and 
30 h"' kpc close pairs, and < Mhost >= 2.2 x 10'- h"' Mq for 
the simulated isolated galaxy sample with Dn > 300 h"' kpc 
and N700 = or 1 . The peaks in the distributions are offset 
by ^0.5 dex. Thus, isolated close pair host halos are ~2 - 3 
times the typical masses of the isolated galaxy (N = 1) sam- 
ple, suggesting that the isolated galaxies are the appropriate 
immediate progenitors of the N = 2 isolated close pairs. In the 



simulation, the subhalos we would identify observationally as 
isolated close pairs were accreted roughly 0-8 Gyr ago, with 
an average of 3.0 Gyr and a wide spread of ct = 1 .5 Gyr For 
comparison, the dynamical timescales for these halos are of 
order 3 Gyr 

As a guide to interpreting this isolated sample, consider 
a simple example. At present, the Milky Way and the An- 
dromeda galaxy, viewed from most projections, would be 
black points on the right side of Fig. |4] They are both iso- 
lated with respect to luminous galaxies. However, as they 
mov e toward one anoth er over the course of the next '-^few 
Gyr dCox & Loebl2007l) . they will move to the left on the plot. 
Eventually their dark matter halos merge and they become an 
N = 2 halo, at which time they will move diagonally toward 
the upper left and "turn green." After perhaps 1-3 Gyr, the 
baryonic galaxies will completely coalesce; the system would 
then move horizontally to the right and "turn black" in the 
plot. Isolated systems like the Milky Way and Andromeda 
are the simplest units of interaction and the best laboratories 
to isolate the effects of interactions from other physical pro- 
cesses. In the next section, we find isolated pairs in the 2dF- 
GRS and use the appropriate control and paired samples to 
isolate the effects of an interaction. As we show in § |7] the 
situation is much more complex in denser environments. 

5. MEASURING TRIGGERED STAR FORMATION 

In § 4, we demonstrate that the amount of triggered star for- 
mation occurring in the universe cannot be measured by com- 
paring typical paired galaxies to typical field galaxies. Dense 
galactic environments are overrepresented in galaxies in pairs. 
However, the effects of an interaction can be identified in- 
dependent from other environmental processes by examining 
isolated pairs and comparing them to isolated galaxies. Al- 
though isolated pairs are not the only environments in which 
triggered star formation is important, they are the simplest to 
identify and study. 

Here, we use 2dF survey pairs to examine triggered star for- 
mation in galaxy pairs. The "full sample" is the full volume- 




FlG. 4. — The distribution of host halo mass as a function of distance to 
the nearest neighbor, Dn, in the lowest-density environments for the objects 
in one realization of our mock galaxy catalogs. The points are the same as 
Fig.|3] but restricted to galaxies with exactly one neighbor within 700 h"' kpc 
and Ay= 1000 km s"' . At right, we plot the host mass histogram of galaxies 
with N700 = 0, or a nearest neighbor distance Dn > 700 h"' kpc. 
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Fig. 5 . — The distribution of host halo masses and separations in low-density environments in our mock catalogs. We repeat Figs.[T]and|2]restricting to galaxies 
with N700 = 1, or exactly one luminous neighbor within 700 h"' kpc. In the left panel, we plot the fraction of galaxies in hosts with a given number of members, 
N, as a function of the distance to the nearest neighbor within 1000 km s"' . At right, we plot histograms. d«/dlog(M), of host halo masses for the control sample 
of N700 < 1 galaxies with Dn > 300 h"' kpc (black, solid) and for < 50 h"' kpc pairs with no other neighbors within 700 h"' kpc (red, dotted) and similarly 
isolated < 30 h"' kpc pairs (blue dashed). The isolated pairs ai'e twice as massive, on average, as the isolated control sample, implying that the control galaxies 
are likely the immediate progenitors of the pairs. 

within 50 h"^ kpc and 1000 km s"' and no others within 700 
h"^ kpc and 1000 km s"'. We also construct the "control" 
sample where N7()o = or 1 and the nearest neighbor is Dn 
> 300 h"' kpc away. 

In the volume-limited sample of 41,239 galaxies, the "full 
sample" includes 22,601 galaxies that are not on the edges 
of the survey volume. There are 1344 galaxies in the "close 
pairs" sample, with spectroscopic companions within 50 h"' 
kpc and 1000 km s"'; 191 of these paired galaxies are in iso- 
lated close (< 50 h~' kpc) pairs and 72 are in isolated closer 
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Fig. 6. — Star-forming properties of 2dFGRS galaxies in pairs. (Top) We 
plot the distribution of rj for the volume-limited sample as a whole (solid, 
black) and for the close pairs sample (dotted, red) with Dn < 5 h~' kpc. 
Large r values of rj con'espond to higher star formation rates I Madg wick et alj 
1200 3°). Relative to the sample as a whole, the pairs exhibit less star formation 
on average, as explained in Fig. [T] by the fact that they preferentially sample 
dense galactic environments. (Bottom) We also restrict to isolated close pairs 
(dotted, red) and a control sample of galaxies with a N700 < 1 and Dn > 300 
h"' kpc distant on the sky. The model shows that the pairs are almost always 
in N = 2 halos (93%) and that 99.5% of the control sample are in N = 1 
halos. Relative to the isolated "control" galaxies, the isolated pairs have more 
galaxies on both ends of the rj distribution. 



limited sample that is far enough away from the edge of the 
survey to sample the environments of galaxies without bias. 
The "close pairs" sample is the subset of this full sample that 
has at least one companion within 50 h"' kpc and 1000 km 
s"'. We focus on the clean measure of the effects of inter- 
actions described in § H by constructing the "isolated close 
pairs" sample of galaxies that have exactly one companion 



(< 30 h" kpc) pairs. The isolated control sample with D^ 
> 300 h"' kpc and N70() = or 1 includes 8564 galaxies. 

We examine the star-forming propert ies of gal a xies u sing 
the spectral parameter ?/, which Madg wick et alJ ('2002) de- 
rive from principal component analysis (PCA) of the 2dF- 
GRS as a whole. Their PCA analysis shows that two-thirds 
of the variance in the 2dFGRS spectra in a volume-limited 
sample are contained in the first two projections, rj is directly 
measured from a combination of these two projections and is 
readily available in the 2dFGRS database. Madgwick et alJ 
(l2003h show that 77 is closely related to both the Ha equiva- 
lent width measured from spectra and to the stellar birthrate 
parameter, b, the present star formation rate of the galaxy di- 
vided by its average past star formation rate. Higher values of 
T? correspond directly t o higher current star formation rates. 
iMadgwick et alJ (l2003h use model spectral energy distribu- 
tions with a wide variety of star formation histories to explore 
this correlation between rj and b. 

In the top of Fig. |6] we plot the distribution of 77 for the 
sample as a whole and all the close pairs. The pairs as a 
whole exhibit, on average, less overall star formation than 
the full volume-limited galaxy sample. The average for the 
full sample and the close pairs sample are < 77 >= 0.07 and 
-0.22, respectively. Fig. [T] demonstrates why. In the model, 
most (56%) galaxies are alone in their dark matter halos. In 
< 50 h~' kpc pairs, however, nearly all (97%) of the galaxies 
are in denser systems, and often much denser systems. Be- 
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Fig. 7. — Nomialized star formation rates of 2dFGRS galaxies in pairs. 
We plot the distribution of b, the ratio of the current star formation rate to 
the average past star formation rate, inferred from r) (see Mad.gwick et al. 
12003 ). We restiict to isolated close pairs {dotted, red) and a control sample of 
galaxies with N700 < 1 and Dn > 300 h"' kpc. A Kolmogorov-Smirnov test 
indicates that the control and close-pair r} distributions differ, with Pk-s = 
6.1 X 10^. Relative to the N = 1 "control" galaxies, the isolated pairs have 
fewer galaxies forming stars at their historical average rate; pairs exhibit an 
excess of galaxies that are very deficient in star formation (5% excess with 
h <0.\) and an excess of galaxies {green shaded region) with rates boosted 
by factors of > 5 (14%) where the high-/) excess, or the shaded region, has 
an average {b) ~ 30. 

cause pairs preferentially reside in extremely dense environ- 
ments, their star formation is suppressed except when they are 
strongly interacting. 

As we demonstrate in §|4l the model predicts that 99.5% of 
the galaxies in the isolated control sample are alone in their 
dark matter halos and that 93% of the isolated pairs are in N 
= 2 halos. Thus, the bottom of Fig.|6]provides us with a nearly 
direct measure of the effects of an interaction, by allowing 
us to compare the true progenitors (the N = 1 "control") of 
the paired (N = 2) systems. A Kolmogorov-Smirnov test indi- 
cates that the control and close-pair -q distributions differ, with 
Pk-s =6.1 x lO"! 

In Fig. [71 we use th e approximate values in Figure 7 of 
iMadgwick et alj (l2003h to convert r] into the ratio of the 
current star formation rate to the average past star formation 
rate, for the isolated control and isolated pairs samples only. 
Relative to the N = 1 "control" galaxies, the isolated pairs 
have fewer galaxies forming stars at their historical average 
rate {b ^ 1). Pairs include an excess of low- and high-star- 
formation-rate galaxies. While 24% of the control sample are 
forming stars with < 0.1, 29% of the < 50 h"' kpc pairs and 
26% of the < 30 h"' kpc are forming stars at these low rates; 
thus, 50 h"' kpc (30 h"' kpc) pairs have a 5% (2%) excess of 
slow (or non-) star formers. 

At the high end of the star-formation rate distribution, 10% 
of the control (N = 1) sample have rates boosted by > a factor 
of 5. In contrast, 50 h"' kpc (30 h"' kpc) pairs are boosted 
by ^7 > 5 24% (30.5%) of the time, for a 14% (20%) excess 
of high star formation rate galaxies. Binning the data with 
fine bins, weighting by b, and limiting the largest (extrapo- 
lated) boost X.ob= 100, the high-/? {b > 5) excess in the pairs, 
the shaded region in Fig. [T] has an average of 31 (34). 
In other words, statistically, this technique shows that for the 
population of pairs with triggered star formation, the average 
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Fig. 8. — Star formation as a function of distance to nearest neighbor in 
sparse environments. {Top) For the 2dFGRS sample with N700 = 1, we plot r] 
as a function of separation to the nearest neighbor; the red line is the averaged 
smoothed by 100 points. The inset displays the average for Dm < 0.13 h"' 
Mpc, highlighting the signal of triggered star formation at Dn < 0.03 h"' 
Mpc. {Bottom) For the N700 = 1 sample in the simulation, we plot the fraction 
of N = 1 {solid, black) and N = 2 {dotted, green) galaxies as a function of Dn , 
as in Fig.|5] 

boost in the star formation rate over the average past rate in 
that galaxy is a factor of ~30. 

The error in the average boost due to triggered star forma- 
tion is dominated by the scatter in the relationship between the 
physical parameter b = SFR/ (SFR) and the measured 77. We 
estimate the error introduced by th i s scatt er using a Monte- 
Carlo simulation. IMadgwick et al.l (l2003h explore the scat- 
ter in the theoretical relationship between b and 77 using syn- 
thetic spectra with star formation histories predicted by semi- 
analytic models. For a typical distribution of star formation 
histories, this scatter is subst antial. We adop t the 1 -a scat- 
ter in the rj-b distribution in IMadgwick et alJ (l2003l) in our 
Monte Carlo simulation, drawing deviations at random from 
the average b for a given (measured) 77. We resample the en- 
tire b distribution for the control, pair, and very close (30 h"' 
kpc) pair rj distributions and remeasure the weighted average 
boost of the excess strong star formers with b> 5 in the pairs. 
Because the scatter of as a function of 77 is much larger in 
the high-/? direction, the 1-cr range of resulting values of the 
average boost of triggered star formation is 42 - 65 for the 
50 h"' kpc pairs and 39 - 61 30 h"' kpc pairs. We conclude 
from this analysis that the large scatter between b and 77 re- 
sults in a wide range of possible star formation boosts result- 
ing from triggered star formation. This type of uncertainty 
would be typical for other popular parameterizations of star 
formation history derived from optical spectra, as well. Many 
optical spectral measures of star formation break down in the 
"bursty" star formation regime that is typical of interacting 
galaxies. 

The excess fraction of starbursting galaxies is enhanced in 
even closer (30 h"' kpc) pairs, while the excess of low star 
formation rate galaxies is not. We explore the dependence 
of starbursting galaxies on separation further by dividing the 
data into bins based on Dn- With 50 h"' kpc bins, only the 
Dn < 50 h"' kpc bin has a statistically significant excess of 
high-/? galaxies compared with the N = 1 control; with 25 h"' 
kpc bins, only the Dn < 50 h"' bins show a significant excess. 
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Fig. 9. — Denser environments: the distribution of host halo mass as a 
function of distance to the nearest neighbor, Dn, but restricted to galaxies 
with 3-8 neighbors within 700 h"' kpc and Ay= 1000 km s"'. The points are 
the same as Fig. [3] 

6. UNDERSTANDING THE TRENDS 

These trends with separation on the sky have been noted 
previously by several authors using rnany different pair 
datasets (Barton et al.'2000':'Lambas et al.'2003V Nikolic et alj 
Ii004; Alonso et al. 2004, 2006; Lin et al. 2007). Thus, the 
fact that we reproduce them here is not surprising. However, 
as Figs. [T] and |5] show, all pair samples suffer from system- 
atic environmental differences as a function of separation on 
the sky that are not directly related to interactions. Typical 
galaxy samples with no environmental controls are dominated 
by cluster galaxies at separations < 200 h"' kpc and field 
galaxies at separations > 200 h~' kpc. Even the most isolated 
samples of pairs are dominated by pairs in N = 2 halos at sepa- 
rations < 200 h"' kpc and by N = 1 halos at separations > 200 
h~' kpc. Here, for the first time, we are uniquely poised to ex- 
plore these effects and construct appropriate control samples 
using models. 

Although we note that isolated galaxies in pairs are the sim- 
plest interactions to study and compare with isolated galaxies, 
they still may be affected by other processes. N = 2 halos are 
often the result of the merger of two N = 1 systems, but they 
can also result from an N = 3 system after two of the three 
baryonic galaxies coalesce. An origin as a N = 3 system is, 
however, much less likely for a given system. Although our 
models are not presently capable of tracking these mergers in 
detail, we note that the density of N = 3 is only 8% of the den- 
sity of N = 1 systems. Thus, the simplest interpretation of the 
result in Fig.|2]is that the 5% excess of low-star-formation sys- 
tems in the N = 2 galaxies results from the fact that galaxies in 
N = 2 halos are more likely than isolated galaxies to have ex- 
perienced previous gas-consuming interactions, and that the 
14% excess of strongly star-forming galaxies in N = 2 halos 
results directly from recent interactions. Future investigation 
of these models will clarify this evolution. 

Supporting this picture, the boosted star formation does not 
extend to all systems with N = 2. In Fig. |8] we plot 77 as a 
function of separation for the isolated (N700 = 1) sample. A 
high mean r\ persists only out to ^-^40-50 h"' kpc. However, 
the model plotted below the points for clarity shows that the N 
= 2 halos dominate well past 100 h"' kpc. Thus, the boosted 



star formation at small separations results from the physical 
interaction and not merely the different mix of halo types at 
smaller separations. A boosted average star formation rate 
does not persist beyond a few tens of h"' kpc, probably be- 
cau se the str ength of the star formation fades as the burst ages 
Csee lBarton et al. 2000: Ba rton Gillespie et al. 2003). and be- 
cause the fraction of galaxies that have had a recent close pass 
is a strong function of Dn. 

The halo substructure model tracks close pericentric pas- 
sages between satellites and the central galaxy in the dark 
matter host halo. These simulations show that for the isolated 
close pair galaxies in true N = 2 systems, 12% of the <50 h"' 
kpc pairs and 16% of the <30 h"' kpc pairs have had a close 
pass with Rperi < 30 h"' kpc within the last Gyr In addition, 
40% of the 50 h"' kpc pairs and 45% of the 30 h"' kpc pairs 
have had a pass with Rperi < 50 h~' kpc within the last 1.2 
Gyr. Thus, the 2dF data are consistent with a range of scenar- 
ios, including at its extremes, (1) the possibility that nearly all 
galaxies with < 30 h"' kpc passes burst by an average factor of 
^30 for 1 Gyr, or (2) the possibility that < half of the galaxies 
with < 50 h"' kpc passes burst this much for ^ 1 .2 Gyr. If star 
formation continues for these long timescales, more widely- 
separated pairs will contain some systems with triggered star 
formation. However, the average star formation rate at wider 
separations is not affected because pairs that have had close 
passes are diluted. The model shows that only 3% of the true 
N = 2 systems with separations 100 < Dn < 200 h kpc have 
had a close (< 30 h"'kpc) pass within the last Gyr 

We will attempt further exploration of these scenar- 
ios in later work. However, preliminary analysis of 
Sloan Digital Sky Survey DR4 NYU Value-Adde d Ga laxy 
Catalog data I York et al ] (120001) : iBlanton et all (HoOSbj); 
lAdelman-McCarthy et alJ (l2006 l) — albeit with quite differ- 
ent star formation measures — suggests that the bursting frac- 
tion is a strong function of the cutoff luminosity of the sam- 
ple studied. The fraction also depends in detail on the star 
formation measure and on the luminosity cutoff used for con- 
structing the volume-limited sample, probably because star 
formation rate is a strong function of galaxy luminosity. Qual- 
itatively, the results remain quite consistent between the 2dF 
and SDSS samples. 

The optical spectra available in the 2dFGRS (and the SDSS) 
do not measure the amount of star formation that is embed- 
ded in dust clouds. Thus, the typical boost of the high-/? ex- 
cess provides an approximate lower limit to the amount of 
trigger ed star formation in a close galaxy-galaxy pass. As 
ICox et al. (2006) demonstrate, high-resolution hydrodynami- 
cal simulations predict a wide range of star formation boosts 
from these close passes. Their models predict SFR boosts 
of 6 - 20 above a galaxy's isolated star formation rate in pro- 
grade encounters during the ^ Gyr between the first close pass 
and the final merger Surprisingly, our estimate of < /> 30 
averaged over the galaxies that have triggered star forma- 
tion exceeds this range. However, the uncertainties associ- 
ated with the mapping between 77 and b are substantial. Thus, 
this 2dF analys is remains broadly consistent with the mod- 
els described in lCox et^D(l2006h and does not yet distinguish 
among them. 

6.1. The effects of survey incompleteness 

The 2dFGRS is not spectroscopically complete; the survey 
is especially deficient in close pairs. We examine the effects of 
incompleteness using the measured completeness of the sec- 
tors in which the galaxies fall, and the 2dF parameters "wsel," 
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a weighted measure of whether the ten nearest galaxies have 
high-quality redshi fts, and "bihrn ," the limiting magnitude o f 
the relevant sector dCoUess etalj|2001l; iNorberg et aLll2002l) . 
Restricting the sample to sectors with known completeness 
> 90% and with wsel > 1.3 and bjlim > 19.35, we find a re- 
duced isolated pairs sample with only 67 pair galaxies and 
2132 control galaxies. As a result, a K-S test no longer reveals 
any statistical significance of the difference in 77 distributions 
(Pk-s = 0.35). Nevertheless, the differences in star-forming 
properties of the isolated control and isolated remain qualita- 
tively unchanged: in the restricted control sample, 12% of the 
galaxies have b> 5;m the restricted isolated pairs, 21% have 
b> 5. Thus, contamination of the isolated control sample ap- 
pears to have no significant effects on our results. 

We further test the effects of redshift survey complete- 
ness by ap plying the s ame analysis to the SPSS pR4 s ample 
dBlanton et al . 2005bi lAdelman-McCarthv et all l2006h . Al- 
though the detailed results are beyond the scope of this paper, 
we are able to use SDSS to find a pure isolated pair sample and 
a pure control sample with SDSS. We construct these sam- 
ples because we are able to exclude all systems with neigh- 
bors that have unmeasured redshifts. The measured excess 
of star formation in isolated pairs compared with the control 
sample, from SDSS depends sensitively on the cutoff lumi- 
nosity of the volume-limited sample. This dependence arises 
because higher-luminosity galaxies have less star formation in 
general. Using Ha equivalent width as a measure of star for- 
mation, the measured excess of pairs with EW(Ha) > 25 A 
ranges from somewhat lower than the 2dF results we present 
here (4%) for a sample volume limited to = -20 +5 log /i to 
19% for a sample volume limited to M,- = -19 + 51og/!. The 
2dF luminosity limit corresponds to galaxies that are midway 
between these numbers. Thus, the 2dF and SDSS analyses 
agree qualitatively, and the fraction of galaxies in close pairs 
with triggered star formation increases when lower luminosity 
galaxies are included. The differences may result from the use 
of very different measures of star formation or from the use of 
galaxies in an intrinsically different luminosity range. How- 
ever, this analysis suggests that "impurity" in the 2dF sample 

— the inclusion of galaxies and pairs that are not truly isolated 

— is not biasing our results significantly. 

7. MORE COMPLEX ENVIRONMENTS 

Isolated pairs of galaxies are the easiest systems to con- 
struct an appropriate control sample for. However, isolated 
pairs are by no means the only pairs that show evidence for 
triggered star formation. Alonso et al.l (120061) show that b still 
increases as pair separation on the sky decreases for close 
pairs in denser environments. 

In Sec. m we describe the effectiveness of isolating the 
lowest-density environments to create an appropriate control 
sample for galaxies in isolated pairs. Unfortunately, the same 
type of technique does not work for higher-density environ- 
ments. One cannot separate galaxies into a density bin and 
expect the non-pairs to be a good control sample for the pairs. 

To illustrate this point, we restrict the sample to interme- 
diate values of our environment statistic; in Fig. |9] we plot 
galaxies with 3 < Nvoo < 8. As expected, isolated N = 1 halos 
are missing from the 50 h"' kpc pairs sample (3.1%) but still 
constitute a significant fraction of the overall sample (31%). 
This result immediately illustrates the problem with break- 
ing sets of galaxies into low, intermediate, and high-density 
samples and separ ating the pairs from these samples (e.g., 
I Alonso et al.ll2006l) . The close pairs have a completely dif- 



ferent environmental mix from the rest of this intermediate- 
density sample. The mix is a strong function of distance to 
the nearest neighbor, and even the pairs reside in a huge range 
of hosts, from h alos with N = 2 galaxies to halos with N > 9. 
However, as Alonso et al.l (1200 6) show, pairs in dense envi- 
ronments do exhibit a rise in star formation rate with a very 
close encounter. 

The ineffectiveness of comparing pairs and "control" galax- 
ies in intermediate- or high-density environments is not 
merely a result of the exact technique we use. Group- 
finding algorithms can be tuned to lower th e incompleteness 
and/or impurity rates on various scales (e.g..lYang et al. 2005t 
Gerke e t al."2005|: |Weinmann et al. l2006l:lBerlind et al.ll2006l: 
Koester et al. 

2007^ 

but these algorithms will always include 
isolated galaxies in denser systems from the control sample at 
a higher rate than they appear in the pairs. Because the distri- 
bution of star formation rates in isolated galaxies has a high- 
star-formation tail, their contamination will always lower the 
measured difference in star formation rate between paired and 
non-paired galaxies in loose groups. One effective approach 
in these environments is to estimate the contamination of N 
= 1 galaxies to the sample of interest and to statistically sub- 
tract its contribution. 

What is the appropriate control sample for pairs in dense 
environments? Even if one could identify the host halo mass 
and occupation of galaxies perfectly, this would still be a dif- 
ficult question to answer. The immediate progenitors of pairs 
in dense environments can range from isolated galaxies that 
have just fallen into the system to galaxies that have been in 
a loose group for along time but are just encountering one an- 
other for the first time. Tracking the evolution of the orbits of 
subhalos in the models will reveal the progenitors of interac- 
tions in crowded environments, but that analysis is beyond the 
scope of the present paper 

8. CONCLUSION 

Here, we use mock catalogs based on simulations of cos- 
mological structure formation in the prevailing ACDM model 
to understand the properties of galaxies in pairs. We examine 
the typical host dark matter halos of galaxies in apparent pairs 
in the simulations and find that: 

1. Our simulations show that galaxies in close pairs are 
preferentially located in cluster and group environ- 
ments. As a result, typical close pairs are not ideal 
for understanding triggered star formation in galaxy- 
galaxy interactions. This result explains why galaxies 
in pairs can appear redder or can appear to have less 
star formation than the typical "field" galaxies. If the 
goal is to isolate the effects of a recent interaction, the 
"field" is not an appropriate control sample with which 
to compare typical galaxies in close pairs. 

2. Using the simulations, we show that close pairs in the 
sparsest environments, with only one neighbor within 
700 h"' kpc, provide the simplest situation for isolating 
the effects of an interaction. Close (< 50 or < 30 h"' 
kpc) isolated pairs are almost exclusively in N = 2 halos 
(93%). Very isolated galaxies with at most one neigh- 
bor within 700 h"' kpc that is at least Dn > 300 h"' kpc 
away are almost exclusively in N = 1 halos (99.5%). 
Thus, isolated field galaxies are the ideal control sam- 
ple of the progenitors of isolated pairs. 
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3. We study isolated pairs and galaxies in a volume- 
limited sample of Mbj < -19 galaxies from the 2dF- 
GRS. Isolated pairs with AD < 50 h"' kpc (< 30 
h~' kpc) show an excess of both strongly star-forming 
galaxies and of non-star-forming galaxies when com- 
pared to N = 1 control sample. While 24% of the con- 
trol sample is forming stars at /?< 0. 1 of its average past 
rate, 29% (26%) of the close pairs lack star formation 
to this extent. In addition, 24% (30.5%) of the pairs are 
forming stars atb>5 times their average past rate while 
only 10% of the control sample has rates this high. The 
rapidly-star-forming excess population in the pairs is 
almost certainly due to the direct effects of the inter- 
action. However, we note that the fraction of galaxies 
undergoing triggered starbursts is a function of the lim- 
iting absolute magnitude of the sample because a much 
smaller fraction of more luminous galaxies are able to 
form stars at all. 

4. For the isolated close galaxy pairs in the 2dFGRS, the 
galaxies with triggered star formation have an average 
star formation boost ofb= SFR/(SFR} - 30. 

5. The measurement of the optical boost of 30 for in- 
teractions in 2dFGRS is approximate, due in large part 
to the large scatter in relationship between star forma- 
tion history and measured optical spectral parameters. 
The measure also excludes embedded star formation 
that is not detectable at optical wavelengths. Thus, it is 
an estimate of the lower limit to the boost in star forma- 
tion rate caused by triggered star formation. This boost 
is higher than predictions derived for boosts from trig- 
gered star formation using hydrodynamic simulations 
(e.g., ICox et al.l 120061) . although the uncertainties are 
too large to rule out specific models at this stage. 



6. By tracking orbits in the substructure model, we show 
that for the isolated close pair galaxies in true N = 2 sys- 
tems, only 12% (16 %) of the 50 (30) h"' kpc pairs have 
had a close pass with Rperi < 30 h"' kpc within the last 
Gyr, and 40% (45%) have had a pass with Rpeii < 50 h"' 
kpc within the last 1.2 Gyr. Thus, the detection of trig- 
gered star formation in 14% (20%) of these L* and sub- 
L* systems in the 2dF data suggest that a large frac- 
tion of the galaxies that experience close passes respond 
with triggered star formation. 
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